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Abstract 

Monitoring  systems  are  proposed  for  the  detection  of  incipient  instability  in  uncer¬ 
tain  nonlinear  systems.  The  work  employs  generic  features  associated  with  the  response 
to  noise  inputs  of  systems  bordering  on  instability.  These  features,  called  “noisy  pre¬ 
cursors”  in  the  work  of  Wiesenfeld,  also  yield  information  on  the  type  of  bifurcation 
that  would  be  associated  with  the  predicted  instability.  The  closed-loop  monitoring 
systems  proposed  in  the  paper  have  several  advantages  over  simple  open-loop  moni¬ 
toring.  The  advantages  include  the  ability  to  influence  the  frequencies  at  which  the 
noisy  precursors  are  observed,  and  the  ability  to  simultaneously  monitor  and  control 
the  system. 


1  Introduction 

In  this  paper,  we  propose  monitoring  systems  for  detection  of  incipient  instability  in  uncertain 
nonlinear  systems.  Our  aim  is  to  develop  automatic  monitoring  systems  that  provide  a 
warning  that  a  system  is  operating  dangerously  close  to  an  instability.  Such  a  warning 
mechanism  can  be  of  great  value  especially  when  no  accurate  system  model  is  available  and 
the  system  is  being  operated  in  a  stressed  condition.  If  an  exact  model  were  available,  then 
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the  stability  boundary  could  be  calculated  off-line,  and  there  would  be  no  need  for  an  on-line 
monitoring  system  for  detecting  incipient  instability. 

Generic-ally,  loss  of  stability  results  in  bifurcation  of  new  steady  states  from  the  nominal 
one  [6],  [17], [19].  The  type  of  bifurcation  that  occurs  depends  on  the  manner  in  which  the 
system  loses  stability.  Since  bifurcations  involve  changes  neglected  in  linearized  models,  it  is 
not  surprising  that  linear  control  design  methods  have  often  been  found  to  be  inadequate  for 
control  law  design  for  stressed  systems.  Bifurcation  control  methods  have  arisen  to  address 
stabilization  of  systems  in  such  situations  [4], 

In  this  paper,  we  contribute  to  this  problem  area  by  considering  the  design  of  monitoring- 
systems  for  detection  of  incipient  instability.  We  also  make  observations  on  the  design  of 
combined  monitoring  and  control  systems  for  physical  systems  susceptible  to  bifurcation  and 
loss  of  stability. 

To  develop  an  on-line  approach  to  the  detection  of  incipient  instability  in  the  absence  of  an 
accurate  system  model,  we  harness  the  effects  of  external  continuously  acting  disturbances. 
The  presence  of  these  disturbances  facilitates  determination  of  measured  signal  features 
associated  with  nearness  to  instability.  The  presence  of  disturbance  inputs,  which  can  occur 
naturally  or  be  injected,  is  crucial.  Without  continuous  disturbances,  a  system  at  equilibrium 
would  remain  at  equilibrium  until  an  instability  occurs,  with  no  possibility  of  an  on-line 
warning  signal.  We  take  the  continuously  acting  disturbances  to  be  white  noise  inputs.  This 
allows  us  to  make  use  of  previous  work  of  Wiesenfeld  and  co-workers  [20,  11,  21,  22,  24,  14,  23]. 
Wiesenfeld  was  interested  in  features  that  can  be  observed  in  the  power  spectrum  of  a 
measured  output  of  a  system  that  is  operating  close  to  an  incipient  bifurcation.  He  focused  on 
nonlinear  systems  operating  along  a  steady  state  limit  cycle.  He  referred  to  the  distinguishing- 
features  in  the  power  spectrum  as  “noisy  precursors”  of  the  bifurcations.  Noisy  precursors 
are  aspects  of  the  power  spectral  density  of  a  measured  output  that  arise  in  the  vicinity  of 
an  instability. 

The  monitoring  systems  we  propose  do  not  employ  or  require  a  system  model.  Rather, 
use  of  the  noisy  precursors  notion  allows  a  nonparametric  approach  based  on  general  features 
of  noise-driven  systems  operating  close  to  instability.  Moreover,  the  monitoring  systems  we 
develop  are  closed-loop,  in  the  sense  that  they  involve  both  sensing  and  actuation.  This 
has  important  advantages  over  the  direct  open-loop  approach  of  simply  monitoring  a  system 
output  and  deciding  if  it  exhibits  features  of  a  noisy  precursor.  Closed-loop  monitoring 
systems  can  enhance  our  confidence  in  deciding  that  system  operation  is  indeed  near  an 
instability  as  well  as  in  determining  the  nature  of  the  instability. 

For  many  engineering  system  models,  the  normal  operating  condition  is  an  equilibrium 
point  rather  than  a  periodic  solution.  Thus,  we  extend  the  theory  of  noisy  precursors  to 
systems  operating  at  an  equilibrium  point.  This  forms  the  basis  for  our  design  of  monitoring 
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systems  for  detection  of  incipient  instability. 

Our  results  apply  to  situations  that  share  the  following  general  characteristics.  A  physical 
system  is  operating  at  or  near  a  nominal  stable  steady  state.  The  system  depends  on  a  set 
of  parameters,  some  of  which  change  slowly  with  time.  Outside  a  certain  range  of  parameter 
values  (the  “design  range”),  model  uncertainty  impedes  reliable  determination  of  system 
stability.  However,  there  are  circumstances  in  which  the  parameters  will  move  outside  of 
the  design  range.  Moreover,  in  these  circumstances  it  is  crucial  that  system  operability  be 
maintained  as  far  as  possible  outside  the  design  range.  The  parameter  changes  may  occur 
due  to  action  of  the  system  operator,  or  may  be  exogenous.  We  refer  to  operation  outside 
the  design  range  as  “stressed  operation.” 

There  are  many  important  examples  in  which  systems  need  to  be  operated  in  off-design, 
stressed  conditions.  The  driving  factors  depend  on  the  application,  but  in  general  they  entail 
a  desire  to  achieve  increased  performance  without  re-design  or  expansion  of  the  system. 
Often,  stressed  operation  leads  to  a  reduced  margin  of  stability.  Thus,  stressed  operation 
can  be  unsafe,  in  that  small  uncertainties  or  disturbances  can  lead  to  loss  of  stability,  i.e., 
to  system  failure.  Examples  include  electric  power  system  voltage  collapse  [18],  chemical 
reactor  runaway  [12],  jet  engine  stall  [13],  aircraft  stall  at  high  angle-of- attack  [5],  and  laser 
system  instability  [8].  For  each  of  these  examples,  precise  models  are  difficult  to  obtain, 
especially  outside  the  design  range. 

The  paper  is  organized  as  follows.  In  Section  2,  we  study  noisy  precursors  for  instability 
for  systems  operating  at  an  equilibrium  point.  In  Section  3,  we  introduce  a  basic  monitoring- 
system  that  facilitates  use  of  precursors  to  detect  incipient  instability.  In  Section  4,  we 
redesign  the  monitoring  system  to  ensure  that  the  bifurcation  occurring  in  the  overall  system 
is  supercritical.  The  monitoring  systems  of  Section  3  and  Section  4  require  full  state  feedback. 
In  Section  5,  we  alleviate  the  full  state  feedback  requirement  for  plants  that  can  be  viewed 
as  singularly  perturbed  (two  time-scale)  systems.  In  Section  6,  we  relax  another  assumption, 
namely  that  the  system  equilibrium  point  is  known.  In  Section  7,  we  give  a  simple  example. 
In  Section  8,  we  collect  our  conclusions. 


2  Noisy  Precursors  for  Nonlinear  System  Instabilities 

In  this  section,  we  extend  the  noisy  precursor  analysis  of  Wiesenfeld  [20]  to  systems  operating 
at  an  equilibrium  point.  Wiesenfeld  considered  systems  driven  by  white  noise  and  operating 
near  a  periodic  steady  state.  He  showed  that  the  power  spectrum  of  a  measured  output 
for  such  a  system  exhibits  sharply  growing  peaks  near  certain  frequencies  as  the  system 
nears  a  bifurcation.  The  particulars  depend  on  the  type  of  bifurcation  that  the  system 
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was  approaching.  He  used  the  results  to  show  that  bifurcating  systems  could  be  used  as 
selective-frequency  amplifiers  [11],  [22],  [23]. 

Consider  a  nonlinear  dynamic  system  ( “the  plant” ) 

x  =  f{x,n)  +  N{t)  (1) 

where  x  G  Rn ,  ^ i  is  a  bifurcation  parameter,  and  N(t)  G  Rn  is  a  zero-mean  vector  white 
Gaussian  noise  process.  Let  the  system  possess  an  equilibrium  point  xq.  For  small  pertur¬ 
bations  and  noise,  the  dynamical  behavior  of  the  system  can  be  described  by  the  linearized 
system  in  the  vicinity  of  the  equilibrium  point  Xq.  The  linearized  system  corresponding  to 
(1)  with  a  small  noise  forcing  N(t.)  is  given  by 

x  =  Df(x0,  fi)x  +  N(t)  (2) 

where  x  :=  x  —  Xq  and  N(t )  G  R"  is  a  vector  white  Gaussian  noise  having  zero  mean.  For 
the  results  of  the  linearized  analysis  to  have  any  bearing  on  the  original  nonlinear  model, 
we  must  assume  that  the  noise  is  of  small  amplitude.  This  assumption  of  small  noise  will  be 
explicated  below,  in  terms  of  smallness  of  correlation  and  cross-correlation  coefficients.  The 
distinct  notation  for  the  system  state  x  and  the  linearized  system  state  x  was  used  here  for 
clarity.  In  the  sequel,  we  will  simply  use  the  notation  x  and  the  meaning  will  be  clear  from 
the  context. 

The  noise  N(t)  can  occur  naturally  or  can  be  injected  using  available  controls.  To  facili¬ 
tate  consideration  of  cases  in  which  the  noise  is  intentionally  injected,  we  write  N(t)  in  the 
general  form 

N{t)  =  Bn(t)  (3) 

where  B  G  Rnxm  and  n(t)  G  Rm  is  a  vector  white  Gaussian  noise.  This  notation  allows 
us  to  easily  consider  cases  in  which  noise  is  injected  in  different  equations  through  available 
actuation  means.  The  noise  vector  N(t)  is  a  white  Gaussian  with  zero  mean  as  long  as  n(t) 
is  a  white  Gaussian  with  zero  mean. 

We  view  the  system  (2)  as  being  in  steady  state  and  driven  only  by  the  noise  process. 
Thus,  we  solve  for  the  evolution  of  the  state  assuming  a  zero  initial  condition.  The  solution 
of  equation  (2)  at  time  t  with  a  zero  initial  condition  is 

x(t)  =  eAt  f  e~AsN{s)cls  (4) 

Jo 

where  A  :=  DF(x o).  For  our  analysis,  we  assume  that  xq  is  an  asymptotically  stable 
equilibrium  point,  i.e.,  all  the  eigenvalues  of  A  have  negative  real  part.  We  can  express  (4) 
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in  terms  of  the  eigenvectors  and  eigenvalues  of  A  (normalized,  and  assumed  distinct): 


x(t) 

zv 


eXkSrkdseXjtrj 


where  r'  and  /'  are  right  and  left  eigenvectors,  respectively,  of  A  corresponding  to  eigenvalue 
A,.  and  where  is  the  Kronecker  delta: 


0  i  ±  j 
1  i  =  j 


(5) 


Thus,  the  i-t.h  component  of  x(t )  is  given  by 


xi(t)  =  Y.^trl  I  lJN(s)e  XjSds 

3= 1 


/0 


(6) 


Since  the  power  spectrum  is  the  Fourier  transform  of  the  autocovariance  function,  we  calcu¬ 
late  the  autocovariance  for  %(t): 


n  n  rt+T  ft 

(xi(t)xi(t  +  r))  =  yyexAex^t+Ar3  k  /  , 

Jo  Jo 


0-\kS2 


j  1  fc=l 

n  n 


•  E  E  (N0{s1)Np(s2))ds1ds2 


o=l  p=l 


Let  the  noise  have  autocorrelation  function 


{.Y;(/!  Y,(/  •  '!)  "<•/>(  r)  (7) 

where  d(-)  is  the  Dirac  delta  function  and  the  iZy  are  constants  for  all  i,j.  Moreover, 
should  be  small  enough  such  that  linearized  analysis  is  valid.  Again,  (7)  is  satisfied  for  the 
linearized  system  as  long  as  n(t)  satisfies 

<>/;{/)//;(/  +  r))  =  4’ijS{r)  (8) 


where  the  'W1;)  are  constants  for  all  i,j.  Then 


(xi{t)Xi{t  +  r)) 


EE 

i=l  k= 1 


eAtexk(t+T)rirk 


rt+r  rt 

Jo  Jo 


e-XjSle~%kS2 
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n  n 

EE  l{lkpvop8(si  -  s2)ds\ds2 

0=1  p=  1 

n  n  pi 

(\,t (\>A<  •  nrirk  /  e~\jS ^—\ks 

1  1  JO 

n  n 

■  E  E  lil>opds  (9) 

0=1  p=l 

Note  that  the  upper  limit  of  integration  changes  from  t  +  r  tot  because  the  impulse  h(si  —  .s2) 
occurs  for  s  1  =  s2,  and  si  <  t  in  the  inner  integral. 

For  a  dynamic  system  that  depends  on  a  single  parameter,  there  are  two  basic  types 
of  bifurcation  from  an  equilibrium  point.  One  is  stationary  bifurcation  in  which  a  new 
equilibrium  emerges  or  the  original  equilibrium  point  suddenly  disappears  at  the  bifurcation. 
The  other  is  Hopf  bifurcation,  where  a  periodic  orbit  emerges  from  the  equilibrium  point 
at  bifurcation.  In  stationary  bifurcation,  a  real  eigenvalue  of  the  linearized  system  becomes 
zero  as  the  parameter  varies.  In  Hopf  bifurcation,  a  complex  conjugate  pair  of  eigenvalues 
crosses  the  imaginary  axis. 

Consider  first  the  Hopf  bifurcation.  Assume  that  a  complex  conjugate  pair  of  eigenvalues 
(denote  them  as  A  =  Ai,  A  =  A2)  close  to  the  imaginary  axis  has  relatively  smaller  negative 
real  part  in  absolute  value  compared  to  other  system  eigenvalues: 


|i2e(A1)|,|i2e(A2)|  <  |i?e(A)| 


(10) 


for  i  =  3, . . . ,  n.  Since  the  integrand  in  (9)  is  the  product  of  decaying  exponentials  (due  to 
the  asymptotic  stability  assumption)  and  bounded,  terms  involving  Ai  and  A2  dominate  (9) 
for  large  t: 


ft  71  71 

<A (/)*,(/  I  r))  eXli2t+T\r])2  <  •A|SEE/!/i'V^ 

J0  7=1  k= 1 


j=ik=i 
n  n 


+  ^,2,+T,(n2)2  /  e-^’YYi’ifads 

j  0  • 17.1 


j=  1  k= 1 


^Ai(t+r)+A2t^ 


rt  n  n 

iri  l  l) 

j  0  • 1  7.  -1 


/.  I-1; jds 


3A2(t+T)+A 


j— 1 k= 1 
n  n 


J0  j=lk=l 


l\vijds 


The  power  spectrum  is  measured  with  the  use  of  a  spectrum  analyzer,  and  most  practical 
spectrum  analyzers  perform  both  an  ensemble  average  and  a  time  average.  Thus,  the  final 
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autocovariance  function  is 


(11) 

(12) 


Cu(t) 


where  ■'  indicates  averaging  over  time  t.  and  we  have  written  C„  in  terms  of  e,  co  >  0  instead 
of  in  terms  of  Ai  =  —  e  +  ju)  and  A2  =  —  e  —  ju.  Also,  H  and  T  are 


Re[{xi(t)xi{t  +  r))  ] 

-  2e'£T  (  \  ,  y  \e~ 

s  • - cos(a>T)  +  I  •  [— 


(ecos(cjr)  —  u  sin(o;r)). 


2(e2  +  u2) 


T 


j= 1  k= 1 

7%  71,  7X  7X 

E  E  +  E  E 

j= 1  &=1  i=l  A:=l 


Finally,  taking  the  Fourier  transform  of  equation  (12)  yields  the  desired  power  spectrum: 


Suiv) 


~  Uv  +  f ) 

"(j/?  +  e)2  + 

T[ _ + f) _ 1 _ 

(e2  +  uj2)((ji]  +  e)2  +  a;2))  (e2  +  co2)((jr]  +  e)2  +  a;2)) 


(13) 


The  magnitude  of  .%;(//!  is  maximum  at  rj  =  uj  and  the  maximum  grows  without  bound  as 
e  — >■  0.  Moreover,  as  the  noise  power  (as  measured  by  the  v^)  increases,  the  magnitude 
of  Sn(rj)  also  increases.  However,  since  H  and  T  affect  .S',;(//)  linearly  and  uniformly  over 
frequency  r 7,  the  shape  of  the  magnitude  Su(ij)  doesn’t  change  with  increasing  noise  power. 
Of  course,  we  have  assumed  that  the  noise  is  of  small  amplitude,  so  we  cannot  actually  allow 
the  Vij  to  increase  without  bound. 

Fig.  1  shows  the  magnitude  of .%,(//)  for  to  =  10,  for  two  values  of  e.  (For  definiteness,  H 
and  T  have  been  set  to  1  in  constructing  Fig.  1.)  Note  the  sharp  peak  around  co  =  10  that 
appears  as  e  — >■  0.  From  this  observation,  we  can  conclude  that  the  power  spectrum  peak 
near  the  bifurcation  is  located  at  u,  and  the  magnitude  of  this  peak  grows  as  e  approaches  to 
zero.  This  property  will  be  used  as  a  precursor  signaling  the  closeness  to  Hopf  bifurcation. 

To  study  the  impact  of  noise  near  a  stationary  bifurcation,  assume  that  a  real  eigenvalue 
close  to  zero  (denote  it  as  A  =  Ai)  and  that  it  has  relatively  smaller  negative  real  part  in 
absolute  value  compared  to  the  other  system  eigenvalues: 


|Ai|  «  |i2e(Ai)| 


(14) 
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for  i  =  2, . . . ,  n.  Due  to  (14),  terms  with  j  =  1  and  k  =  1  dominate  the  expression  (9)  for 
large  t.  so  that 

n  n 

j= 1 k= 1 

Taking  the  time  average,  we  get  the  autocovariance  function 

Cu(r)  :=  (xi{t)Xi{t  +  t))*  (15) 

=  EE^teKh1)2^-  (16) 

7  I  k= 1  “  e 

Fourier  transformation  of  (16)  gives  the  desired  power  spectrum: 

n  n  i 

W  =  (IT) 

This  equation  shows  that  the  magnitude  of  the  power  spectrum  peak  grows  as  e  approaches 
to  zero  and  the  location  of  this  peak  is  rj  =  0.  Fig.  2  shows  the  magnitude  of  Sri  (  rj)  (17).  (For 
definiteness,  the  coefficient  in  square  brackets  in  Eq.  (17)  has  been  set  to  1  in  constructing 
Fig.  2.)  Note  the  sharp  growing  peak  around  lo  =  0  as  e  — >■  0. 

3  Monitoring  System  for  Detecting  Incipient 
Stationary  and  Hopf  Bifurcation 

As  shown  in  the  foregoing  section,  we  can  expect  to  observe  a  growing  peak  in  the  power 
spectrum  of  a  measured  output  of  a  nonlinear  system  with  white  Gaussian  noise  input  as  the 
system  approaches  a  bifurcation.  In  the  case  of  Hopf  bifurcation,  the  location  of  the  power 
spectrum  peak  coincides  with  the  imaginary  axis  crossing  frequency  of  the  critical  eigenval¬ 
ues.  In  the  case  of  stationary  bifurcation  the  power  spectrum  peak  occurs  at  zero  frequency. 
In  this  section,  we  use  these  observations  to  develop  a  monitoring  system  for  proximity  to 
bifurcation.  Since  noisy  precursors  associated  with  stationary  bifurcation  involve  a  grow¬ 
ing  peak  in  the  power  spectrum  at  zero  frequency,  these  are  difficult  to  resolve.  Hence,  we 
first  propose  a  closed-loop  monitoring  system  that  addresses  this  problem  by  transforming 
a  stationary  bifurcation  into  a  Hopf  bifurcation.  That  is,  the  original  plant  augmented  with 
the  monitoring  system  undergoes  a  Hopf  bifurcation.  The  critical  frequency  of  the  Hopf 
bifurcation  is  set  by  the  monitoring  system  itself.  After  introducing  the  monitoring  system 
and  studying  its  use  in  monitoring  for  stationary  bifurcation,  we  study  its  use  in  monitoring 
a  system  for  proximity  to  a  Hopf  bifurcation. 


{Xi{t)Xi{t  +  t)) 


DXi(2t+r) 


(rift 

Jo 
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3.1  Generating  a  Hopf  Bifurcation  from  a  Stationary  Bifurcation 

Suppose  the  plant  of  interest  is  susceptible  to  loss  of  stability  through  a  stationary  bifur¬ 
cation.  Since  Hopf  bifurcation  is  easier  to  detect  than  stationary  bifurcation  through  noisy 
precursors,  we  introduce  a  monitoring  system  that  replaces  the  stationary  bifurcation  with 
a  Hopf  bifurcation  of  tunable  frequency. 

In  the  absence  of  noise,  let  the  plant  obey  the  dynamics 

&  =  f{x,fi)  (18) 

Results  we  obtain  for  this  model  will  have  immediate  implications  for  precursor-based  mon¬ 
itoring  of  the  system  with  noise  effects  included.  Suppose  the  following  assumptions  hold: 

(51)  The  origin  is  an  equilibrium  point  of  system  (18)  for  all  values  of  //,. 

(52)  System  (18)  undergoes  stationary  bifurcation  from  the  origin  at  //  =  jic  (i.e.,  there  is 
a  simple  eigenvalue  A (y)  of  Df(0,y)  such  that  for  some  value  //  =  //.r. ,  A(//r:)  =  0  and 

dMp-c)  Q) 

(53)  All  other  eigenvalues  of  Df( 0,  jic)  are  in  the  open  left  half  complex  plane. 

We  introduce  the  following  augmented  system  (plant  plus  monitoring  system)  correspond¬ 
ing  to  (18): 


Xi  =  f(x,  fi)  —  cyi 

Vi  =  cXi  (19) 

Here,  y  G  Rn ,  c  G  R  and  /  =  1,2,...,  n.  Eq.  (19)  will  later  be  viewed  as  a  basic  monitoring- 
system  whose  use  facilitates  detection  of  either  stationary  or  Hopf  bifurcation.  Note  that 
the  state  vector  consists  of  the  original  physical  system  states  x  augmented  with  the  states 
y  of  the  monitoring  system. 

Proposition  1  Under  assumptions  (SI)- (S3),  the  augmented  system  (19)  undergoes  a  Hopf 
bifurcation  from  the  origin  at  y  =  yc.  Moreover,  if  for  any  value  of  y  the  origin  of  the  original 
system  (18)  is  asymptotically  stable  (resp.  unstable),  then  the  origin  is  asymptotically  stable 
(resp.  unstable)  for  the  augmented  system  (19). 
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Proof:  Denote  by  A  the  Jacobian  matrix  of  system  (18)  at  the  origin.  Clearly,  the  origin 
(0,0)  in  R2n  is  an  equilibrium  point  of  the  augmented  system  (19).  The  Jacobian  matrix  of 
the  augmented  system  (19)  at  the  origin  is 


A  -cl 
cl  0 


(20) 


Let  a  be  any  eigenvalue  of  A  and  r  the  corresponding  right  eigenvector.  Also,  denote  by  A 
any  eigenvalue  of  J  and  the  associated  right  eigenvector  by  v  =  [tq  v2]T .  Then, 


Aui  =  Avi  —  cv2 

(21) 

Xv2  =  cv  1 

(22) 

We  seek  a  solution  for  which  v\  =  r.  From  (22),  we  have 

c 

,2  = 

(23) 

Substituting  (23)  into  (21)  and  using  r  ^  0,  we  get 

A2  -  aA  +  c2  =  0 

(24) 

Thus,  any  eigenvalue  a  of  A  has  corresponding  to  it  two  eigenvalues  of  which  are  the 
solutions  of  the  quadratic  equation  above: 

a  zb  \/ a2  —  4c2 

A  =  V  - —  (25) 

2 

Thus,  the  eigenvalues  of  the  Jacobian  matrix  of  the  augmented  system  (19)  are 

a,:  ±  \/a?  —  4c2 

A2?:— 1,2?:  = - ^ -  i  =  1,  2, . . . ,  n 

where  a,;,  i  =  1, . . . ,  n  are  the  eigenvalues  of  A.  Let  the  eigenvalue  of  A  that  becomes  0  at 
criticality  be  aq.  At  //  =  jic,  the  eigenvalues  of  the  augmented  system  associated  with  ai  are 
(using  (24))  are  a  pair  of  pure  imaginary  eigenvalues  at  c 

Ai,A  2  =  ±cj  (27) 

Note  that  the  pair  of  pure  imaginary  eigenvalues  (27)  depends  on  c. 
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For  a  Hopf  bifurcation  to  occur,  the  transversality  condition  should  be  satisfied.  That 
is,  the  eigenvalues  crossing  the  imaginary  axis  should  do  so  with  nonzero  speed.  From  the 
quadratic  equation  (24)  and  using  the  fact  that  oc\  =  0  at  /i  =  jic,  we  have 

dRe{ Ai)  _  1  cloii 

dji  2  dfi 

Since  cq  =  0  and  ^  ^  0  at  //  =  jic  from  assumption  (S2),  (28)  implies  dRed^  =  ^  0 

(i.e. ,  the  transversality  condition  holds  for  system  (19)).  Therefore,  the  augmented  system 
(19)  undergoes  a  Hopf  bifurcation  from  the  origin  at  /i  =  jic. 

The  last  step  in  the  proof  consists  in  showing  that  all  other  eigenvalues  of  the  matrix 
.J  are  in  the  open  left  half  complex  plane.  Any  pair  of  eigenvalues  of  J  can  be  obtained 
from  (26).  For  a  real  eigenvalue  of  A ,  it  is  clear  from  (26)  that  the  corresponding  pair  of 
eigenvalues  of  J  have  negative  real  part  if  cq  <  0  since  a.,  <  Re{  \J a 2  —  4c2  }.  For  a  complex 
conjugate  pair  of  eigenvalues  of  A  (denoted  7,7),  we  have  the  following  two  equations: 


A2  —  7  A  +  e2  =  0  (29) 

A2  — 7A  +  C2  =  0  (30) 

Multiply  (29)  and  (30)  to  get  the  following  fourth  order  equation: 

A4  -  (7  +  T)A3  +  (2c2  +  77) A2  -  c2(7  +  7)A  +  c4  =  0  (31) 

Denoting  7  =  a  +  bj  and  7  =  a  —  bj,  equation  (31)  simplifies  to 

A4  -  2a A3  +  (2c2  +  a2  +  b2)\2  -  2ac2A  +  c4  =  0  (32) 


Applying  the  Routh-Hurwitz  criterion  [7]  to  (32),  we  obtain  the  Routh  array 


•s4  1  2c2  +  a2  +  b2  c4 

S'3  —2  a  —2  ac2  0 

s 2  3  c2  +  a 2  +  b 2  c4  0 

1  —2  ac2  (2c2 -\-a62-\-b2) 

5  3c2+a2+62  U  U 

s°  c4  0  0 


From  assumption  (S3),  a  <  0.  This  guarantees  that  all  the  entries  in  the  first  column  of  the 
Routh  array  are  positive.  Therefore,  all  eigenvalues  of  the  Jacobian  matrix  of  the  augmented 
system  have  negative  real  part. 
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From  the  foregoing  discussion,  it  is  also  clear  that  if  any  eigenvalue  of  A  has  positive 
real  part,  then  the  corresponding  eigenvalues  of  J  also  have  positive  real  part.  This  proves 
that  if  the  origin  is  unstable  for  the  plant,  then  it  is  also  unstable  for  the  augmented  system.  ■ 

Note  that  since  the  value  c  in  equation  (19)  is  adjustable,  we  can  control  the  crossing- 
frequency  of  the  complex  conjugate  pair  of  eigenvalues  of  the  augmented  system.  Thus, 
for  detecting  stationary  bifurcation,  we  only  need  to  monitor  a  frequency  band  around  the 
chosen  value  of  c.  It  is  also  possible  to  slowly  vary  c  in  a  controlled  fashion,  giving  added 
confidence  in  our  assessment  that  an  instability  is  imminent. 

There  are  some  other  advantages  of  our  monitoring  system.  The  augmented  system  (19) 
has  the  same  critical  parameter  value  (yc)  as  the  original  system.  This  is  actually  not  a 
luxury  but  a  necessity  for  the  system  to  be  practically  useful.  In  addition,  the  final  part 
of  the  proof  shows  that  augmenting  the  states  yi  and  applying  the  feedbacks  cyi  to  the 
original  system  does  not  change  the  local  stability  of  the  system.  Moreover,  to  apply  the 
monitoring  system,  we  do  not  need  knowledge  of  the  original  system.  However,  there  are 
some  restrictions  and  further  considerations  in  applying  the  suggested  monitoring  system  to 
general  physical  systems.  We  will  discuss  these  in  subsequent  sections. 

We  assumed  above  that  the  stationary  bifurcation  was  such  that  the  transversality  con¬ 
dition  is  satisfied.  This  means  that  the  eigenvalue  that  vanishes  at  criticality  must  cross 
into  the  right  half  of  the  complex  plane  with  nonzero  speed  as  the  parameter  is  varied.  In  a 
saddle  node  bifurcation,  however,  the  nominal  equilibrium  disappears  at  criticality,  so  that 
the  transversality  condition  does  not  hold.  For  a  saddle-node  bifurcation,  the  augmented 
system  of  this  section  results  in  a  degenerate  Hopf  bifurcation.  The  possible  bifurcation 
diagrams  for  degenerate  Hopf  bifurcation  are  more  complex  than  for  Hopf  bifurcation  [16], 
[10].  However,  for  the  purpose  of  detecting  incipient  instability,  the  details  of  the  ensuing  bi¬ 
furcation  are  not  important.  These  details  become  important  when  we  consider  the  system’s 
post-bifurcation  behavior. 

3.2  Detecting  Hopf  Bifurcation  using  Monitoring  System 

In  this  section,  we  consider  the  effect  of  the  monitoring  system  of  the  proceeding  section  on 
a  system  that  undergoes  Hopf  bifurcation  instead  of  stationary  bifurcation.  Consider  again 
the  system  (18),  repeated  here  for  convenience: 

x  =  f(x,y)  (33) 
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(HI)  The  origin  is  an  equilibrium  point  of  (33)  for  all  values  of  //,. 

(H2)  System  (33)  undergoes  a  Hopf  bifurcation  from  the  origin  at  //  =  jic. 

(H3)  All  other  eigenvalues  of  Df( 0,  //r.)  are  in  the  open  left  half  complex  plane. 

As  in  the  foregoing  section,  let  the  augmented  system  (plant  plus  monitoring  system)  be 

Xi  =  Mx,n)-qji 

Vi  =  cXi  (34) 

where  x  G  Rn ,  y  G  Rn ,  cG  1?,  and  /  =  1,2,...,  n. 

Proposition  2  Under  the  assumptions  (H1)-(H3),  the  augmented  system  (34)  undergoes  a 
codimension  two  bifurcation  at  p  =  yc,  in  which  two  complex  conjugate  pairs  of  eigenvalues 
cross  the  imaginary  axis.  Moreover,  for  any  value  of  fi  if  the  origin  of  the  original  system  is 
asymptotically  stable  (resp.  unstable),  then  the  origin  is  asymptotically  stable  (resp.  unstable) 
for  the  augmented  system. 

Proof:  First,  we  show  that  the  augmented  system  has  two  pairs  of  pure  imaginary  eigenval¬ 
ues  at  the  origin  for  //,  =  fic,  and  that  these  eigenvalues  satisfy  the  transversality  condition. 

From  assumption  (H2),  the  Jacobian  matrix  of  the  original  system  at  the  origin  has  a 
pair  of  pure  imaginary  conjugate  eigenvalues  (denote  them  by  ju) ,  —juf)  for  //,  =  pc.  From  the 
proof  of  Proposition  1,  it  is  clear  that  each  of  these  eigenvalues  results  in  a  pair  of  eigenvalues 
for  the  augmented  system  which  are  the  solutions  of  the  following  equations: 

A2  —  ju)\  +  c2  =  0  (35) 

A2  +  jco  A  +  c2  =  0  (36) 

By  multiplying  the  equations  above,  we  get  a  fourth  order  equation  the  solutions  of  which 
are  eigenvalues  of  augmented  system: 

A4  +  (2c2  +w2)A2  +  c4  =  0  (37) 

The  four  solutions  of  the  equation  above  are  given  by 
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Note  that  2c2  +  uj2  >  V 4 c2co2  +  ad  for  all  c.  x  G  R.  Therefore,  the  Jacobian  matrix  of 
the  augmented  system  has  two  pairs  of  pure  imaginary  conjugate  eigenvalues  at  the  critical 
parameter  value. 

To  check  the  transversality  condition,  consider  the  eigenvalues  for  //  near  ^ ic .  Near  g  =  gc. 
we  have  the  following  fourth  order  equation  the  solutions  of  which  result  from  the  pair  of 
complex  conjugate  eigenvalues  (a  ±  uj)  of  the  original  system  (see  (31)): 

A4  -  2aA3  +  (2c2  +  «2  +  u2)\2  -  2 c2a  +  c4  =  0  (39) 


Since  g  is  close  to  /ic,  by  continuity  it  follows  that  this  equation  has  two  pairs  of  complex 
conjugate  eigenvalues  as  its  solutions  for  //  near  gc.  Denote  these  as  e±fj,g±hj.  The  next 
relationship  is  now  easily  demonstrated: 

4 

-2a  =  ^A  i  =  e  +  g  (40) 

i—  1 

—2c2 a  =  ^2  =  2 g{e2  +  / 2)  +  2.e{g2  +  h 2)  (41) 

i,j,k=l 


where  A,-.  Ar  A/,,  are  roots  of  (39).  Taking  the  derivative  of  both  sides  of  the  equation  above 
with  respect  to  //  and  evaluating  at  //  =  /i.c  (e  =  g  =  0,  j3  =  0,y2  =  c  —  a)  gives 


h 2 


i  de 

dji 


de  dg 
dg  dg 

f2dl 

dg 


r)  do. 
dg 
2  da 
dg 


(42) 


Next  we  solve  Eq.  (42)  for  Also,  h  and  /  are  not  0  at  the  critical  point  from  (38) 

and  /  /  /;  at  the  critical  point  if  c  7^  0.  These  conditions  guarantee  that  if  ^  ^  0,  then 

%  ^  ^  °- 

The  last  step  in  the  proof  consists  in  showing  that  all  other  eigenvalues  of  the  Jacobian 
matrix  J  of  (34)  lie  in  the  open  left  half  complex  plane.  Note  that  we  have  the  same  form 
of  matrix  J  as  in  the  proof  of  Proposition  1: 


A  -cl 
cl  0 


(43) 


where  all  noncritical  eigenvalues  of  A  have  negative  real  part.  We  can  use  the  same  procedure 
as  in  Proposition  1  to  prove  that  if  all  noncritical  eigenvalues  of  A  have  negative  real  part, 
then  all  corresponding  eigenvalues  of  J  have  negative  real  part. 
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It  is  also  clear  that  if  any  eigenvalue  of  A  has  positive  real  part,  then  the  corresponding 
eigenvalues  of  J  also  have  positive  real  part.  This  implies  that  if  the  original  system  is 
unstable,  then  the  augmented  system  is  also  unstable.  ■ 

Since  two  pairs  of  eigenvalues  of  the  augmented  system  cross  the  imaginary  axis  at  the 
critical  parameter  value,  we  can  expect  to  see  two  peaks  in  the  power  spectrum  as  the  system 
nears  the  bifurcation  point.  From  (38),  we  see  that  the  values  of  the  pairs  of  imaginary 
eigenvalues  at  criticality  depend  on  c.  Hence,  we  can  change  the  location  of  the  power 
spectrum  peaks  by  changing  c.  Moreover,  we  can  predict  the  exact  locations  of  the  peaks  if 
the  pair  of  eigenvalues  crossing  the  imaginary  axis  in  the  original  system  is  known. 

Because  two  pairs  of  eigenvalues  cross  the  imaginary  axis  for  the  augmented  system,  the 
augmented  system  undergoes  a  codimension  two  bifurcation.  The  nature  of  the  bifurcation 
behavior  depends  strongly  on  The  possible  bifurcation  diagrams  for  the  associated 

degenerate  Hopf  bifurcation  can  be  found  in  [16],  [10].  However,  as  was  the  case  for  de¬ 
tecting  incipient  stationary  bifurcation,  the  details  of  the  degenerate  Hopf  bifurcation  are 
not  important.  They  become  important  when  we  consider  the  system’s  post-bifurcation 
behavior. 


4  Stabilization  of  Bifurcated  Limit  Cycle 
in  the  Augmented  System 

In  this  section,  we  suppose  that  the  plant  is  subject  to  loss  of  stability  through  a  stationary 
bifurcation.  In  these  circumstances,  the  monitoring  system  proposed  above  results  in  a  Hopf 
bifurcation  in  the  augmented  system.  Besides  being  able  to  predict  that  a  bifurcation  is 
about  to  take  place,  it  would  be  useful  if  the  monitoring  system  could  also  ensure  stability 
of  the  bifurcated  solution.  That  is,  a  system  that  can  perform  both  monitoring  and  control 
functions  is  desirable.  The  purpose  of  this  section  is  to  illustrate  how  the  monitoring  system 
we  have  proposed  can  be  modified  to  serve  in  both  capacities.  Liberal  use  is  made  of  the 
bifurcation  formulas  and  associated  results  summarized  in  Appendix  B. 

4.1  Stability  of  the  Bifurcated  Limit  Cycle 
of  the  Augmented  System 

First,  we  consider  the  relationship  between  the  stability  of  bifurcated  equilibrium  points  of 
the  original  system  and  stability  of  the  bifurcated  limit  cycle  of  the  augmented  system.  If 
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stability  of  the  bifurcated  equilibria  of  the  original  system  implies  stability  of  the  bifurcated 
periodic  solution  of  the  augmented  system,  then  bifurcation  control  design  need  only  be 
performed  for  the  original  system.  We  proceed  to  show,  however,  that  the  stability  properties 
of  the  bifurcation  persist  in  the  case  of  scalar  systems,  but  not  generally  for  systems  of 
dimension  two  or  higher. 

Let  the  plant  be  given  by 

X  =  f(x1fi)  (44) 

where  x  G  R"  is  the  state  vector  and  y  €  R  is  the  bifurcation  parameter,  and  the  noise  input 
is  neglected  for  the  purposes  of  this  section.  Suppose  that  at  the  critical  parameter  value 
//  =  fic,  the  Jacobian  matrix  of  (44)  evaluated  at  the  equilibrium  point  xo  =  0  has  one  zero 
eigenvalue. 

For  simplicity,  we  first  consider  the  one-dimensional  case  (n  =  1),  i.e.,  suppose  the 
state  vector  of  (44)  is  a  single  variable.  Also,  suppose  system  (44)  undergoes  a  pitchfork 
bifurcation.  It  is  easy  to  see  that  the  left  (l)  and  right  (r)  eigenvectors  corresponding  to 
the  simple  zero  eigenvalue  at  criticality  can  be  taken  as  any  nonzero  constants.  Set  r  =  1 
and  l  =  1,  so  that  r  and  l  satisfy  the  normalization  Ir  =  1.  Lemma  1  (see  Appendix  B) 
then  applies  directly,  allowing  calculation  of  the  associated  bifurcation  stability  coefficients 
pi  and  p2-  As  discussed  in  Appendix  B,  the  pitchfork  bifurcation  is  supercritical  (giving- 
stable  bifurcated  equilibria)  if  P\  =  0  and  82  <  0.  Since  system  (44)  is  assumed  to  undergo 
a  pitchfork  bifurcation,  Pi  =  0: 


Pi  =  lQ(r ,  r)  =  ^(0)  =  0  (45) 

Thus,  Q(r,r)  vanishes.  Since  Q(r,r )  =  0,  we  have  x2  =  0  (using  the  notation  of  Appendix 
B).  Thus,  p2  becomes 

2  d3  f 

P2  =  2lC0(r,  r,  r)  =  — ^-(0)  (46) 

The  augmented  system  corresponding  to  the  plant  (44)  is 

x  =  f  (x,  n)  —  cy 

y  =  cx  (47) 

where  x,  y,  //  G  R.  We  have  shown  that  the  augmented  system  undergoes  a  Hopf  bifurcation  if 
the  original  system  undergoes  a  pitchfork  bifurcation.  To  check  the  stability  of  the  bifurcated 
periodic  solution  of  (47),  we  only  have  to  check  the  sign  of  the  Hopf  bifurcation  stability 
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coefficient  /T  (102).  At  criticality,  the  Jacobian  matrix  of  (47)  is 


L0 


0 

c 


(48) 


The  matrix  L0  has  an  eigenvalue  cj  with  corresponding  right  eigenvector  r  =  [  1  —  j  j  and 
left  eigenvector  l  =  |  [  1  j  j .  Eigenvalue  —  cj  has  right  eigenvector  f  and  left  eigenvector  l. 
Note  that  higher  order  terms  only  come  from  f(x ,  //)  and  they  are  not  a  function  of  y.  From 
this  observation,  we  have 


1  h  [ x  v] 

Q5|Q3 

tor-* 

3 

0  " 

X 

\ 

Q{{x,y),{xiy))  = 

[0 

0 

. y . 

\ 

0 

/ 

ifA0)*2  1 


(49) 


From  equation  (45),  equation  (49)  implies  that  Q{r,  r )  and  Q(r,  r)  both  vanish.  Therefore, 
the  solutions  of  (100)  and  (101)  are  a  =  0  and  6  =  0.  Now,  /T  of  (47)  becomes 

4  3  1  f)2  f 

f32  =  2Re{-lC(r,r,f)}  =  2--^(0)  (50) 

since  higher  order  terms  only  come  from  f(x,[i)  and  none  depend  on  y. 

Note  that  the  sign  of  (46)  agrees  with  the  sign  of  (50).  The  next  proposition  therefore 
follows. 


Proposition  3  Suppose  the  system  (44)  is  of  first  order,  i.e.,  n=l.  If  the  plant  (44)  un¬ 
dergoes  a  supercritical  pitchfork  bifurcation  (respectively  a  subcritical  pitchfork  bifurcation) , 
then  the  transformed  system  (47)  undergoes  a  supercritical  Hopf  bifurcation  (respectively  a 
subcritical  Hopf  bifurcation). 

Next,  we  consider  the  case  n  >  2,  that  is,  the  case  in  which  the  dimension  of  the  plant 
is  at  least  2.  We  show  using  an  example  that  the  monitoring  system  proposed  above  does 
not  necessarily  preserve  the  stability  character  of  the  bifurcation  in  the  plant.  That  is,  a 
supercritical  pitchfork  bifurcation  (resp.  a  subcritical  pitchfork  bifurcation)  in  the  plant 
need  not  result  in  a  supercritical  Hopf  bifurcation  (resp.  a  subcritical  Hopf  bifurcation)  in 
the  augmented  system. 
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Consider  the  example 


X\  =  —fix  1  —  x\  +  XiX2 

X-j  =  -x2-\-kx\  (51) 

where  //  G  i?  is  a  bifurcation  parameter  and  A:  G  R  is  a  constant.  It  is  easy  to  see  that  the 
origin  is  an  equilibrium  point  for  all  parameter  values  /i  and  that  a  pitchfork  bifurcation 
occurs  for  fi  =  0.  Moreover,  a  simple  calculation  shows  that  fi2  for  this  pitchfork  bifurcation 
is  —1.  This  implies  that  system  (51)  undergoes  a  supercritical  pitchfork  bifurcation  at  y  =  0. 
The  augmented  system  corresponding  to  (51)  is 

x\  =  —fix  i  —  x\  +  a;  1^2  —  q/i 

2/i  =  c.xi 

x2  =  —x2  +  kx\  -  cy2 

y2  =  cx2  (52) 

As  discussed  in  Appendix  B,  typically  a  Hopf  bifurcation’s  stability  is  determined  by  a  single 
bifurcation  stability  coefficient  ,d2  (this  differs  from  the  fi2  coefficient  in  the  study  of  pitchfork 
bifurcations).  The  Hopf  bifurcation  is  supercritical  if  the  coefficient  j32  is  negative,  and  it  is 
subcritical  if  the  coefficient  is  positive.  We  now  calculate  fi2  for  the  Hopf  bifurcation  that 
occurs  in  the  augmented  system  (52).  To  facilitate  application  of  the  formulas  in  Appendix 
B,  denote  the  state  vector  of  (52)  as  ,2  =  ( zi,  Z2,  23,  Za)T  where  z\  :=  X\,  z2  :=  yi,  23  :=  x%, 
and  24  :=  y2.  The  Jacobian  matrix  of  (52)  evaluated  at  the  origin  at  criticality  is 

'  0  -c  0  0 

c  0  0  0 

0  0  -1  -c 

0  0  c  0 


(53) 


One  eigenvalue  of  this  matrix  is  cj,  and  it  has  corresponding  right  eigenvector  r  =  [  1  —j  0  0  j 
and  left  eigenvector  1  =  |  1  j  0  0  j .  The  conjugate  eigenvalue  —cj  has  right  eigenvector 

f  and  left  eigenvector  l.  The  Taylor  series  expansion  of  the  right  side  of  (52)  has  the  following- 
quadratic  and  cubic  terms: 


Q(z,z) 


'  Z1Z3  ' 

X\X2 

0 

0 

kz\ 

hr2 

IvJy  2 

0 

0 
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C(z,z,z) 


(54) 


/O- Jb 

0  _  0 

0  “  0 

0  J  0 

Therefore,  we  have 

'  0  ' 

Q(r,r)  =  Q(r,r)=  ° 

.  0  . 

'  -1  ' 

C(r,r,r)  =  (55) 

0 

Solving  Eqs.  (100)  and  (101)  of  Appendix  B,  we  obtain 

a  =  [  0  0  0  -f  ]T 

&  =  [0  0  rfirc  Asf  («) 

Substituting  these  values  in  (102)  of  Appendix  B,  we  find  that  /T  for  system  (52)  is 


*  =  -I  +  IT9? 

Note  that  for  sufficiently  large  k ,  /T  is  positive.  For  such  values  of  k,  the  augmented  system 
(52)  therefore  undergoes  a  subcritical  Hopf  bifurcation  even  though  the  plant  undergoes  a 
supercritical  pitchfork  bifurcation.  Thus,  for  n  >  2,  the  monitoring  system  as  presented 
above  does  not  necessarily  preserve  the  stability  of  bifurcated  solutions.  We  now  proceed  to 
modify  the  design  to  address  this  deficiency. 


4.2  Redesign  for  Combined  Monitoring 
and  Bifurcation  Stabilization 

In  this  section,  we  modify  our  monitoring  system  such  that  the  bifurcated  limit  cycle  occur¬ 
ring  in  the  augmented  system  is  guaranteed  stable  regardless  of  the  stability  of  the  pitchfork 
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bifurcation  occurring  in  the  plant.  That  is,  we  ensure  that  the  Hopf  bifurcation  in  the 
augmented  system  is  supercritical,  and  that  this  holds  regardless  of  whether  the  pitchfork 
bifurcation  in  the  original  system  is  supercritical  or  subcritical.  The  modification  that  we 
introduce  in  the  monitoring  system  involves  the  addition  of  a  nonlinear  term  with  a  gain 
parameter  that  can  be  tuned  to  ensure  the  desired  result. 

Let  the  plant  obey  the  dynamics 


x=f{x,n),  (58) 

which  we  assume  undergoes  a  pitchfork  bifurcation  from  the  origin  for  n  =  jic.  Here,  x  G  R" 
and  //  G  R  is  the  bifurcation  parameter.  Denote  by  rs  and  ls  the  right  and  left  eigenvectors, 
respectively,  corresponding  to  the  simple  zero  eigenvalue  of  the  system  linearization  at  //  = 
He.  Take  the  first  component  of  rs  to  be  1  and  impose  the  normalization  lsrs  =  1  (following 
the  procedure  in  Appendix  B). 

Now  consider  the  following  redesign  of  the  augmented  system: 

Xi  =  fi{x,  n)  —  ci/i 
Vi  =  cXi  -  mxlvi 

Here,  c  and  m  are  real  constants.  At  criticality,  the  system  (59)  has  Jacobian  matrix 

T_  \  A  -cl 
'  ~  [  cl  0 

at  the  origin,  written  in  terms  of  A,  the  Jacobian  matrix  of  (58).  Employing  Proposition  1, 
it  is  easy  to  show  that  the  augmented  system  (59)  undergoes  Hopf  bifurcation  and  that  the 
matrix  J  has  eigenvalues  ±ci.  Moreover,  the  right  and  left  eigenvectors  of  J  corresponding 
to  the  eigenvalue  cj  are  given  by 

r  =  [  rs  —jrs  ]T 

l  =  \[h  jls  ]  (61) 

Also,  the  right  and  left  eigenvectors  corresponding  to  the  eigenvalue  —  cj  are  given  by  f  and 
I,  respectively. 

The  stability  of  the  bifurcated  periodic  solution  of  the  augmented  system  is  determined 
by  the  sign  of  the  bifurcation  stability  coefficient  /A  (Eq.  (102),  Appendix  B): 

$2  =  2Re{2lQo(r,a)  +  lQ0(r,b)  +  A'i /  .  i  r  I  [  (62) 


(60) 


(59) 
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Since  there  are  no  quadratic  terms  in  in  the  augmented  system  (59),  /32  simplifies  to 

=  2Re{2lsQf(r,  a)  +  lsQf{r ,  b )  +  ic(r,  r,  r )}  (63) 

where  Q/(-,-)  denotes  the  quadratic  terms  in  the  Taylor  expansion  of  f(x,nc).  Moreover, 
we  can  simplify 


lC{r,  r,  r)  =  l8Cf{r ,  r,  r)  +  J-lsCy{r,  r,  r) 

m  n 

=  ls(7J(r,r,f)  +  -^K(rs1)2  (64) 

D  i=i 

where  Cf(x,x,x)  denotes  the  cubic  terms  in  the  Taylor  expansion  of  f(x,nc),  and  r\  and  lls 
denote  the  i-tli  component  of  rs  and  ls,  respectively.  Since  the  first  component  of  r]  is  1  and 
lsrs  =  1,  (64)  reduces  to 

771 

lC(r,  r,  r)  =  lsCf{r ,  r,  r)  -  —  (65) 

Hence,  /32  becomes 

3  3777, 

02  =  2Re{2lsQf(r,  a)  +  lsQf{r,  b )  +  -IJ  >(/•.  r,  r)}  -  —  (66) 

By  choosing  m  positive  and  sufficiently  large,  we  can  ensure  that  ,d2  will  be  negative.  This 
will  imply  that  the  Hopf  bifurcation  occurring  in  the  augmented  system  (59)  is  supercritical. 

Here,  we  have  suggested  only  one  of  many  possible  designs  that  render  the  Hopf  bifurcated 
supercritical.  The  method  is  robust,  since  the  efficacy  of  the  design  does  not  depend  on  the 
details  of  the  plant  model.  For  any  given  plant,  a  sufficiently  large  feedback  gain  m  will 
result  in  supercriticality  of  the  Hopf  bifurcation.  Note  that  we  added  a  nonlinear  term  only 
to  the  dynamics  of  the  augmented  states  y,  not  to  those  of  the  physical  system  states  Xj.  We 
therefore  have  considerable  freedom  in  choosing  the  nonlinear  feedback  gain  m. 

5  Reduced  Order  Monitoring  System 

The  closed-loop  monitoring  systems  introduced  in  the  preceding  two  sections  entail  the  use 
of  full  state  feedback.  In  this  section,  we  alleviate  this  requirement  for  plant  models  that 
can  be  viewed  as  singularly  perturbed  (or  two  time-scale)  systems.  We  design  a  monitoring- 
system  in  which  only  the  slow  states  are  fed  back  to  the  controls. 
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Consider  a  plant  given  by  a  singularly  perturbed  system  of  the  form 

x  =  f{x,  z,  n,  e) 

ez  =  g{x,z,(i,e)  (67) 

where  x  G  Rn,  z  G  Rm ,  //.  e  E  R  and  e  is  small  but  positive.  The  reduced  system  is  obtained 
by  formally  setting  e  =  0  in  (67),  giving 

oc  =  /(./'.  z,  fi.  0) 

0  =  g{x,  z,  fi,  0)  (68) 

Let  mo  =  (0,zo)  be  an  equilibrium  point  of  the  reduced  system.  Also,  assume 
(SP1)  m.0  =  (0,  Zq)  is  an  equilibrium  point  of  (68)  for  all  values  of  //,. 

(SP2)  f,g ,  are  Cr(r  >  5)  in  .t,  2.  //,,  e  in  a  neighborhood  of  (m,o,  0,  0). 

(SP3)  No  eigenvalue  of  Z)2(/(0,  2o,  0,  0)  has  zero  real  part. 

(SP4)  The  reduced  system  undergoes  a  stationary  bifurcation  at  m.0  for  the  critical  parameter 
value  f.i  =  //r.. 

Let  the  augmented  system  (plant  plus  monitoring  system)  corresponding  to  (67)  be 

x  =  f{x,  z,  fi.  e)  -  cy 
ij  =  cx 

ez  =  g(x,  z,  n,  e)  (69) 

Proposition  4  Let  (SP1)-(SP4)  above  hold.  Then  there  is  an  e0  >  0  and  for  each  e  G  [0,  e0] 
the  augmented  system  (69)  undergoes  a  Hopf  bifurcation  at  an  equilibrium  mjc’f  near  m-o  for 
a  critical  parameter  value  pec  near  pc. 

Proof:  By  virtue  of  Theorem  in  [1]  on  persistence  of  Hopf  bifurcation  under  singular  pertur¬ 
bation,  we  need  to  verify  two  conditions.  The  first  is  that  the  reduced  system  corresponding 
to  (69)  undergoes  a  Hopf  bifurcation  at  (0,  0,  zq)  at  the  critical  parameter  value  //  =  pc.  The 
second  condition  is  that  the  Jacobian  matrix  of  g  with  respect  to  the  fast  variables  2  does 
not  possess  any  eigenvalues  with  zero  real  part.  The  reduced  system 

x  =  f{x,  z,  fi,  0)  -  cy 
y  =  cx 

0  =  g{x,z,n,  0)  (70) 
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Since  the  original  reduced  system  (68)  undergoes  a  stationary  bifurcation,  we  can  apply 
Proposition  1  to  (70)  to  show  that  the  reduced  augmented  system  (70)  undergoes  Hopf  bi¬ 
furcation  at  the  critical  parameter  value  //  =  jic.  The  result  now  follows  from  [1].B 

Proposition  4  is  useful  because  it  implies  that  we  only  have  to  augment  and  feed  back  slow 
states  in  a  two-time  scale  system  to  transform  stationary  bifurcation  into  Hopf  bifurcation. 

6  Monitoring  System  for  Nonzero  Equilibrium  Point 

Although  the  results  of  the  preceding  sections  do  not  depend  on  availability  of  an  accurate 
model  of  the  plant,  they  do  require  knowledge  of  the  nominal  equilibrium.  In  this  section,  we 
alleviate  this  requirement  through  a  re-design  of  the  monitoring  system.  Not  surprisingly, 
the  increased  generality  comes  with  some  cost,  mainly  in  the  simplicity  of  the  observed 
instability  precursor. 

The  requirement  of  a  known  equilibrium  is  embodied  in  assumption  (SI),  which  states 
that  the  nominal  equilibrium  point  of  the  plant  is  fixed  at  the  origin  for  all  parameter  values. 
Through  a  simple  parameter-dependent  coordinate  change,  it  is  clear  that  the  results  of  the 
preceding  sections  still  apply  under  the  milder  assumption  that  the  equilibrium  is  a  known 
function  of  the  parameter. 

Assumption  (SI)  was  invoked  so  that  the  equilibrium  point  of  the  plant  is  not  changed 
upon  applying  state  feedback.  A  standard  control  technique  for  exactly  preserving  an  equi¬ 
librium  despite  model  uncertainty  involves  the  use  of  washout  filters  [15].  Our  revised  designs 
in  this  section  entail  adjoining  a  washout  filter  to  the  previous  monitoring  system  designs. 

Denote  by  Xq[/i.)  the  nominal  equilibrium  point  of  system  (18).  We  now  allow  the  equi¬ 
librium  to  depend  in  some  unknown  fashion  on  the  parameter  /i. 

The  following  re-designed  augmented  system  involves  two  sets  of  additional  variables: 
the  vector  y ,  which  appears  in  the  original  design  (19);  and  the  vector  z,  the  washout  filter 
states: 


Xi  =  fiix./i)  ci/i 
Vi  =  cXi  +  azi 

A  =  Vi  (71) 

for  i  =  1,2,...,  n,  where  a,  c  G  R. 

Proposition  5  Assume  the  original  system  (18)  satisfies  (S2)  and  (S3)  at  an  equilibrium 
point  x0 {n)  not  necessarily  at  the  origin.  Then  the  augmented  system  (71)  undergoes  a 
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codimension  2  bifurcation  at  fi  =  fic.  At  criticality,  the  linearization  of  (71)  possesses  one 
simple  zero  eigenvalue  and  a  pair  of  pure  imaginary  eigenvalues. 

Proof:  The  equilibrium  point  of  the  augmented  system  (71)  is  (x0,  0,  z 0),  where  z0  is  solution 
of  cXj+ciZj  =  0.  Note  that  new  augmented  system  keeps  x0  as  a  component  of  this  equilibrium 
point.  The  Jacobian  matrix  of  (71)  evaluated  at  this  equilibrium  point  is 

A  —cl  0 

J  =  cl  0  al 

0  /  0 

where  A  is  the  Jacobian  matrix  of  the  original  system  evaluated  at  x0.  Let  a  be  any  eigen¬ 
value  of  A  and  r  corresponding  eigenvector.  Also,  assume  A  is  an  eigenvalue  of  J  with 
eigenvector  v  =  [vf  vff  vf]T.  Then 


If  we  choose  a  <  0,  then  J  has  eigenvalues  0,  ±Vc2  —  aj  which  correspond  to  the  zero 
eigenvalue  of  the  original  system  criticality. 
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Next,  we  check  the  transversality  condition.  Equation  (76)  which  corresponds  to  crossing 
simple  real  eigenvalue  of  original  system  has  one  real  and  a  pair  of  complex  conjugate  as 
its  solution  near  the  critical  point.  Denote  8  as  the  real  eigenvalue  and  [3  ±  yj  as  the 
pair  of  complex  conjugate  eigenvalue.  Solving  this  notation  in  (76)  and  separating  real  and 
imaginary  parts,  we  obtain 


5  -\-  2  *  j3  =  —a 
8{/32  +  y2)  =  act 

Differentiating  these  equations  with  respect  to  //,  gives 

d<)  ^  d(3 

dji  djt 

^{fi  +  ^)  +  AS  +  2%S  = 

djt  djt  dji 


da 

djt. 


da 

djt. 


(78) 


(79) 


At  the  critical  parameter  value  //,  =  jic,  8  =  0,  f3  =  0,  and  y 2  =  c2  —  a.  Thus,  at  //  =  /y. 


dS  d/3  da 

- h  2—  = - 

dji  dji  dji 

dS  2  da 

y (c  “ (,)  = 

From  equation  (81),  4^  ^  0.  Solving  these  equations  for  gives 


(80) 

(81) 


d/3  1  c  da 

dji  2  c2  —  a  dji 


(82) 


which  is  nonzero  if  7^  0  and  a  <  0. 

As  was  the  case  with  Proposition  1,  the  final  step  in  the  proof  consists  of  showing  that 
all  other  eigenvalues  of  the  matrix  J  are  in  the  open  left  half  complex  plane  (C_).  There  are 
three  eigenvalues  of  J  which  correspond  to  one  negative  real  value  eigenvalue  of  A  and  these 
eigenvalues  are  solutions  of  the  (76).  By  using  the  Routh-Hurwitz  criterion,  we  can  show 
that  solutions  of  equation  (76)  in  C_  if  the  corresponding  real  eigenvalue  of  A  is  in  C_.  For 
the  complex  conjugate  pair  of  eigenvalues  of  A  (y,y),  we  have  following  two  equations 


A3  —  yA2  +  (c2  —  a)  A  +  ay  =  0 
A3  —  yA2  +  (c2  —  a)A  +  ay  =  0 


(83) 

(84) 
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Multiply  (83)  and  (84)  to  get  the  sixth  order  equation 

A6  -  (7  +  y)A5  +  (2(c2  -  a)  +  77)A4  +  (7  +  y)(2a  -  c2)A3 

+  ((c2  —  a)2  —  2a77)A2  +  (c2  —  a)a(  7  +  7)  A  +  a2  77  =  0  (85) 

By  applying  the  Routh-Hurwitz  criterion  to  (85),  we  can  show  that  all  solutions  of  (85)  are 
in  the  open  left  half  complex  plane  if  Re( 7)  is  negative  (details  are  in  Appendix  A).  ■ 

We  have  proved  that  the  new  augmented  system  (71)  with  nominal  equilibrium  not 
necessarily  at  the  origin  replaces  a  stationary  bifurcation  with  a  codimension  two  bifurcation. 
Note  that  the  design  gives  the  same  critical  parameter  value  for  the  plant  and  the  augmented 
system.  In  addition,  the  crossing  eigenvalues  at  critical  point  are  located  0  and  ±\/c2  —  aj. 
Also,  note  that  an  original  simple  zero  eigenvalue  persists  under  the  augmentation.  Thus,  the 
monitoring  system’s  effectiveness  has  to  do  with  its  introduction  of  a  purely  imaginary  pair 
of  eigenvalues  at  criticality  in  addition  to  the  zero  eigenvalue.  Near  bifurcation,  we  expect 
power  spectrum  peaks  to  be  located  at  0  and  (?  —  a.  By  varying  c  and  a  (both  tunable 
parameters),  we  can  tune  the  location  of  the  peak  at  \/c2  —  a  as  desired.  This  flexibility 
increases  our  assurance  that  the  power  spectrum  peak  is  caused  by  closeness  to  instability 
rather  than  by  other  factors  (such  as  noise).  However,  the  new  augmented  system  (71)  also 
comes  with  some  disadvantages  compared  to  the  system  in  Proposition  1.  In  Proposition  1, 
we  transform  a  stationary  bifurcation  into  a  Hopf  bifurcation.  In  other  words,  the  system 
has  a  limit  cycle  as  its  solution  instead  of  a  new  equilibrium  point  near  of  bifurcation.  In 
comparison  to  the  previous  augmented  system  design  (19),  the  new  augmented  system  (71) 
shows  more  complicated  bifurcation  behavior  [9] .  The  system  is  no  longer  guaranteed  to  have 
a  periodic  orbit  as  a  solution  near  bifurcation.  Either  a  periodic  orbit  or  a  new  equilibrium 
point  could  result  at  bifurcation.  The  bifurcation  diagram  depends  strongly  on  the  vector 
field  f(x,n).  However,  it  may  be  possible  that  augmented  system  has  desired  bifurcation 
diagram  by  introducing  some  nonlinear  terms  into  augmented  states.  Of  course,  to  do  that 
we  have  detail  knowledge  on  f(x,n).  Details  on  codimension  two  bifurcations  can  be  found 
in  [9].  However,  for  the  purpose  of  monitoring,  it  is  enough  to  have  a  discernible  power 
spectrum  peak  when  the  system  approaches  instability. 

The  next  proposition  asserts  that  the  new  augmented  system  also  works  for  singularly 
perturbed  systems  using  fewer  states  for  feedback.  The  only  difference  from  the  previous 
results  on  singularly  perturbed  systems  is  that  we  no  longer  require  (SP1)  of  Section  5. 
Using  the  same  notation  as  in  Section  5,  we  have  the  following  proposition. 

Proposition  6  Let  (SP2)-(SP4)  of  Section  5  hold  for  the  system  (18).  Then  there  is  an 
e0  >  0  and  for  each  e0  £  [0,  e0]  the  following  extended  system  undergoes  a  codimension  two 
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(one  real  and  a  pair  of  complex  eigenvalues  crossing)  bifurcation  at  an  equilibrium  m!(c'e  near 
?77-o  for  a  critical  parameter  value  yec  near  yc: 


Xi 

=  fi(x,  /i.<)  ciy 

Vi 

=  cxj,  +  awi 

m 

=  Vi 

ez 

A 

d 

II 

(86) 

where  i  =  1,2, ,  n. 

Proof:  Follows  directly  from  Proposition  5  and  Theorem  7  of  [1].  ■ 

7  An  Example 

Consider  again  the  simple  system  (51)  which  undergoes  a  pitchfork  bifurcation.  For  conve¬ 
nience,  we  rewrite  the  equations  for  the  plant  (including  a  noise  term)  augmented  with  a 
monitoring  system: 

X\  =  —fix  i  —  x\  +  X\X2  —  cyi  +  N(t) 

Vi  =  cx  i 

x 2  =  -x2  +  5xj  -  cy2 

y2  =  cx2  (87) 

Here,  N(t)  is  a  white  Gaussian  noise.  The  system  (87)  undergoes  a  Hopf  bifurcation  at 
//,  =  0,  which  is  the  parameter  value  where  a  pitchfork  bifurcation  occurs  for  the  original 
system  (51).  The  origin  loses  stability  as  //,  is  decreased  through  //,  =  0. 

The  simulation  results  in  this  section  were  obtained  by  the  MATLAB  Simulink  package. 
Figure  3  shows  the  location  of  the  power  spectrum  peak  in  frequency  as  the  parameter  c 
is  (quasistatically)  changed.  The  simulations  were  done  for  a  parameter  value  of  y  =  0.1, 
which  is  before  the  origin  loses  stability.  Note  from  the  figure  that  the  location  of  the  power 
spectrum  peak  obtained  from  simulation  (shown  as  an  asterisk  in  Figure  3)  agrees  well  with 
the  predicted  location  (straight  line  in  Figure  3). 

8  Conclusions 

We  have  proposed  closed-loop  monitoring  systems  for  detection  of  incipient  instability  in 
uncertain  nonlinear  plants.  These  systems  make  use  of  characteristics  of  the  power  spectrum 
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of  a  measured  output  in  the  vicinity  of  an  instability.  By  employing  closed-loop  designs,  we 
are  able  to  more  reliably  monitor  for  incipient  instability  through  on-line  tuning  of  control 
parameters.  Two  time-scale  models  were  used  to  reduce  the  number  of  measurements  fed 
back  in  the  closed-loop  monitoring  system.  We  have  studied  the  impact  of  the  monitoring 
systems  on  stability  of  bifurcations  that  occur  when  stability  is  lost,  and  proposed  design 
modifications  to  ensure  stability  of  bifurcated  solutions  in  a  robust  fashion. 

A  Routh-Hurwitz  Calculation  for 
Proposition  5 

Using  the  same  notation  as  in  Proposition  5,  and  letting  7  =  a  +  ;jco  and  7  =  a  —  juj,  Eq. 
(85)  becomes 


A6  -  2<jA5  +  (2(c2  -  a)  +  a2  +  uj2) A4  +  2<r(2a  -  c2)A3  + 
{{c2  —  a)2  —  2a{a2  +  uj2))\2  +  2{c2  —  a)aa\  +  ci2{a2  +  uj2)  =  0 


Applying  the  Routh-Hurwitz  criterion  to  the  equation  above,  we  obtain  the  Routh  array 


1 

-2a 

c  T  a2  T  us2 

2ac(a—  (<72+cj2)) 
c+a2+u)2 

A 

— 2ac3a(c+a2+L02) 
a(a2+uj2+c)-(a+c)2 

ci2  (a2  +  u)2) 


2  (c  —  a)  +  a  ( c 

2a(2a  —  c ) 

c 2  —  ac  —  2  a(a2  +  to2) 

2aac(c—a.+a2  +lo2) 

C+CT2+(jJ2 

a2(a2  +  cj2) 

0 

0 


a)2  —  2  a(a2  +  uj2)  a2(a2  +  uj2) 


2  aa(c-a)  0 

a2  (a2  +  cj2)  0 

0  0 

0  0 

0  0 

0  0 


where  A  =  U^+^iCc+^+^-Ca+c)2^2^2)  If  a  <  0  and  a  <  0,  then 

a—{a2+w2)  ’ 

column  of  this  array  are  positive.  Therefore,  under  this  condition  all 
negative  real  part. 


all  entries  in  the  first 
solutions  of  (85)  have 


B  Stability  of  Bifurcated  Solutions 

In  this  appendix,  we  recall  formulas  for  determining  the  local  stability  of  bifurcated  solutions. 
Details  can  be  found  in  Abed  and  Fu  [2], [3]. 
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Consider  a  one-parameter  family  of  nonlinear  autonomous  systems 

x  =  f(x,fi)  (88) 

where  x  G  Rn  is  the  vector  state  and  //,  is  a  real-valued  parameter.  Let  f(x,  //)  be  sufficiently 
smooth  in  x  and  //  and  let  :ro,yu  be  the  nominal  equilibrium  point  of  the  system  as  a  function 
of  the  parameter  //,. 

First,  we  consider  the  case  of  stationary  bifurcation.  For  simplicity,  we  take  the  critical 
parameter  value  to  be  jic  =  0  in  the  statement  of  the  next  hypothesis. 

(S)  The  Jacobian  matrix  of  system  (88)  at  the  equilibrium  x0^  has  a  simple  zero  eigenvalue 
Ai (n)  with  A'^O)  ^  0,  and  the  remaining  eigenvalues  lie  in  the  open  left  half  of  the 
complex  plane  for  //  =  0. 

The  Stationary  Bifurcation  Theorem  asserts  that  hypothesis  (S)  implies  a  stationary 
bifurcation  from  Xo4Ji  at  fi  =  0  for  (88).  A  new  equilibrium  branch  bifurcates  from  x^4J  at 
H  =  0.  The  theorem  states  that  near  the  point  (:ro,o,0)  of  the  (n  +  l)-dimensional  (xu u)- 
space,  there  exists  a  locally  unique  curve  of  critical  points  {x(e),  yu(e)),  distinct  from  x0 and 
passing  through  (.To,o,0),  such  that  for  all  sufficiently  small  |e|,  x(e)  is  an  equilibrium  point 
of  (88)  when  //,  =  ju(e).  (Here,  e  is  an  auxiliary  small  parameter.) 

The  series  expansions  of  x(e),/u(e)  can  be  written  as 

yU(e)  =  /LllC  +  yU2f2  +  •  •  •  (89) 

x(e)  =  x0^  +  X\t  +  x2e2  +  •  •  •  (90) 

If  /i i  ^  0,  the  system  undergoes  a  transcritical  bifurcation  from  x0^  at  //  =  0.  That  is, 

there  is  a  second  equilibrium  point  besides  Xq^  for  both  positive  and  negative  values  of  // 

with  |//|  small.  If  ji\  =  0  and  //2  ^  0,  the  system  undergoes  a  pitchfork  bifurcation  for  |//| 

sufficiently  small.  That  is,  there  are  two  new  equilibrium  points  existing  simultaneously, 
either  for  positive  or  for  negative  values  of  //  with  //  small.  The  new  equilibrium  points 
have  an  eigenvalue  /3(e)  which  vanishes  at  //  =  0.  The  series  expansion  /3(e)  is  given  by 


/3(e)  —  pie  +  /32e2  +  •  •  • 

(91) 

We  have  the  exchange  of  stability  formula: 

Pi  =  — A7(0) 

(92) 

Moreover,  in  case  pi  =0,  p2  is  given  by 

P2  =  —  2^2A7(0) 

(93) 
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Eqs.  (92), (93)  are  not  explicit  formulas  for  pi  and  p2.  Explicit  formulas  are  given  Lemma  1. 

The  bifurcation  stability  coefficients  P\  and  p2  can  be  obtained  using  eigenvector  com¬ 
putations  and  series  expansion  of  the  vector  field.  System  (88)  can  be  written  in  the  series 
form 


x  =  Lpx  +  Qp\x,  x)  +  Cn(xy x.x)  I - 

=  Lqx  T  q iL\x  qp  L %x  T  •  •  • 

+Q0(x,x)  +  nQ \{x,x)  - 

I  C\)(.r.x.x)  I -  (94) 

where  x  =  x  —  xo,o,  L^,  L\,  L2  are  n  x  n  matrices,  Qn(x,  x),  Qo(x,x),  Qi(x,  x)  are  vector¬ 
valued  quadratic  forms  generated  by  symmetric  bilinear  forms,  and  C^lx,  x,  x),  Co(x,  x,  x) 
are  vector-valued  cubic  forms  generated  by  symmetric  trilinear  forms. 

By  assumption,  the  Jacobian  matrix  L0  has  a  simple  zero  eigenvalue  with  the  remain¬ 
ing  eigenvalues  stable.  Denote  by  l  and  r  the  left  (row)  and  right  (column)  eigenvectors, 
respectively,  of  the  matrix  L0  associated  with  the  simple  zero  eigenvalue,  where  the  first 
component  of  r  is  set  to  1  and  the  left  eigenvector  l  is  chosen  such  that  Ir  =  1.  (Setting 
the  first  component  of  r  to  1  sometimes  requires  a  re-ordering  of  the  state  variables.)  The 
following  well  known  fact  is  used  in  the  statement  of  the  next  lemma: 

A'(0)  =  lLxr  (95) 

The  two  lemmas  that  follow  give  stability  criteria  for  the  bifurcated  equilibria  of  system 
(94).  The  first  addresses  to  pitchfork  bifurcation,  while  the  second  addresses  transcritical 
bifurcation. 

Lemma  1  Let  hypothesis  (S)  hold.  Then 


Pi  =  lQo{r,r) 

(96) 

Also,  if  pi  =  0,  then 

p2  =  2l{2Q0{r,  x2)  +  Co(r ,  r,  r)} 

(97) 

where  X2  solves 

L0x  2  =  ~Qo(r,r) 

(98) 

The  bifurcated  equilibrium  points  of  (94)  near  x Q  for  p.  near  0 
pi  =  0  and  p2  <  0,  but  they  are  unstable  if  pi  =  0  and  /J2  >  0. 

are  asymptotically  stable  if 

33 


Lemma  2  Let  hypothesis  (S)  hold,  and  suppose  that  (3\  ^  0  (this  can  be  checked  using  (96)). 
Then  the  bifurcated  solution  is  asymptotically  stable  on  one  side  of  p  =  0  and  is  unstable  on 
the  other.  For  any  given  value  of  p  near  0,  the  stability  of  the  bifurcated  solution  is  opposite 
that  of  the  nominal  equilibrium. 

Now  consider  system  (88)  under  for  the  following  hypothesis,  which  implies  occurrence  of 
Hopf  bifurcation.  Again,  for  simplicity  the  critical  value  of  the  parameter  is  taken  as  jic  =  0. 

(H)  The  Jacobian  matrix  of  system  (88)  at  the  equilibrium  x0^=0  has  a  pair  of  pure  imag¬ 
inary  eigenvalues  Ai(0)  =  jujc  and  Ai(0)  =  —  jujc  with  coc  ^  0,  the  transversality 
condition  dReW0)}  ^  0  is  satisfied,  and  all  the  remaining  eigenvalues  lie  in  the  open  left 
half  complex  plane. 

Under  these  conditions,  the  Hopf  Bifurcation  Theorem  asserts  the  existence  of  a  one- 
parameter  family  pe,  0  <  e  <  eo  of  nonconstant  periodic  solutions  of  system  (88)  emerging 
from  x  =  xo )At  at  the  parameter  value  0  for  sufficiently  small  \p\.  Exactly  one  of  the  charac¬ 
teristic  exponents  of  p,  governs  the  asymptotic  stability  and  is  given  by  a  real,  smooth  and 
even  function 

/3(e)  =  /32e2  +  /tf^e4  +  •  •  •  (99) 

Specifically,  pe  is  orbitally  stable  if  /3(e)  <  0  but  is  unstable  if  /3(e)  >  0.  Generically  the  local 
stability  of  the  bifurcated  periodic  solution  pe  is  decided  by  the  sign  of  the  coefficient  /32. 
It  happens  that  the  sign  of  /32  also  determines  the  stability  of  the  critical  equilibrium  point 
x0 An  algorithm  for  computing  the  stability  coefficient  /32  follows. 

Stepl  Express  (88)  in  the  Taylor  series  form  (94).  Let  r  be  the  right  eigenvector  of  L0 
corresponding  to  eigenvalue  jcoc  with  the  first  component  of  r  set  to  1.  Let  l  be  the  left 
eigenvector  of  L0  corresponding  to  the  eigenvalue  jcoc,  normalized  such  that  Ir  =  1. 

Step2  Solve  the  equations 

L0a  =  ~^Q0(r,f) 

(2 jujcI  -  L0)b  =  i Q0(r,r ) 

for  a  and  b. 

Step  3  The  stability  coefficient  /32  is  given  by 

/32  =  2Re{2lQ0(r,  a)  +  IQo(r.  b )  +  ^/G(r,  r,  r)}  (102) 


(100) 

(101) 
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